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Summary 

The unusually low frequency oscillation in the wake of an airfoil, studied in 
[1], is explored experimentally as well as computationally for a NACA0012 airfoil 
with a "glaze ice accretion" at the leading edge. Experimentally, flow oscilla- 
tions are observed at low frequencies that correspond to a Strouhal number of about 
0.02. This occurs in the angle of attack range of 8° to 9°, near the onset of 
static stall for this airfoil. With a Navier-Stokes computation, "limit-cycle" 
oscillations in the flow and in the aerodynamic forces are also observed at low 
Strouhal numbers. However, the occurrence of the oscillation is found to depend 
on the turbulence model in use as well as the Reynolds number. 

Nomenclature 
c airfoil chord 

Ci lift coefficient 
Cd drag coefficient 

C m moment coefficient about 0.5c 

M Mach number 

f s shedding or oscillation frequency 

R c chord Reynolds number 

St s Strouhal number, f ? c sina/U 

u'(f) one-dimensional, longitudinal velocity spectrum 
U free-stream mean velocity 

a angle of attack 


1 . Introduction 


An unusually low frequency oscillation in the flow over a LRN ( 1 ) - 1 007 airfoil 
was studied experimentally as well as computationally in [1], over the R c range 
of 15x10-3 to 300x10-3. The phenomenon involved a naturally occurring periodic 
switching between stall and unstall around the onset of the static stall condition. 
The frequency was considered low as the corresponding Strouhal number, based on 
the cross-stream length scale and the free stream velocity, was only about 0.02, an 
order of magnitude lower than that experienced in the commonly observed bluff-body 
shedding (see also [2]). Even though the mechanism of frequency selection has 
remained unresolved, the measurements established that the phenomenon was indeed 
aerodynamic in origin. Any connection to a standing acoustic wave, structural res- 
onance or unsteadiness associated with the tunnel was ruled out. The flow oscilla- 
tions imparted very large unsteady forces to the airfoil. Thus, investigations 
leading to a clearer understanding of the phenomenon are well in order from the 
point of view of applications involving unsteady aerodynamic loads, e.g., in stall 
flutter of blades and wings. 

In [1], unsteady oscillations were also observed computationally that had 
striking similarities with the experimental results. In an independent study [3], 
symptoms of the phenomenon were also encountered in a Navi er-Stokes computation 
for the flow over a NACA0012 airfoil with "glaze ice accretion." A glaze ice 
accretion develops during flight conditions just below the freezing temperature. 
These shapes are typified by large "horns" (see the simulated shape in Figure 1(b)). 
This led to a wind tunnel experiment with an airfoil model having the same cross- 
sectional shape as used in [3]. Lift characteristics showed the onset of static 
stall around a » 7° with this airfoil. Slightly above this a, wake oscillations 
were observed resulting in an unambiguous peak in the velocity spectra, although 
not as sharp as in the "LRN" airfoil case, at a Strouhal number of about 0.02! 

The objective of the present paper is to summarize these results. Salient 
features of the phenomenon observed with the "LRN" airfoil are reviewed first. 

The experimental results obtained with the "iced" airfoil are then presented, fol- 
lowed by a discussion of limited computational results. As it will be observed, 
even though the computation seems to capture the essence of the phenomenon, the 
results vary widely depending on the turbulence model in use and the Reynolds 
number . 
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2. Experimental Procedure 

The experiments were carried out in a low speed wind tunnel which has been 
described in detail elsewhere [1,4]. Figure 1(a) is a photograph of the test sec- 
tion as set up during a somewhat similar experiment [4]. The test section is 76 cm 
(wide) by 51 cm (high). The maximum speed is about 16 ms -1 , however, some of the 
data were obtained using a smaller drive motor yielding a maximum speed of about 
10 ms -1 . The free-stream turbulence intensity is less than 0.1 percent, but could 
be raised to about 0.4 percent by inserting a 30-mesh screen about 21 cm upstream 
of the airfoil support. Two-dimensional airfoil models were used in the experi- 
ment. Data from the LRN ( 1 ) - 1 007 and a NACA0012 airfoil with simulated "glaze ice 
shape" at the leading edge are presented. Both are of nominal chord of 12.7 cm 
and of aspect ratio of about 6. These two airfoils, whose cross sections are 
shown in Figure Kb), are referred to simply as the "LRN" and the "iced" airfoils, 
respectively. The photograph in Figure 1(a) shows the "LRN" airfoil mounted in 
the tunnel . 

Standard hot-wire anemometry was used for velocity measurements. A Nicolet 
660B analyzer was used for spectrum analysis. There was provision for acoustic 
excitation of the flow, and an automated angular positioning device was used dur- 
ing the measurement of the lift and the drag [4]. The airfoils were supported at 
both ends at mid-chord. The coordinate origin is at the airfoil mid-chord and the 
axis of the test section, x, y and z denote streamwise, vertical and spanwise 
coordinates, respectively. The angle of incidence a is measured with respect to 
the (dashed) chord lines as shown in Figure Kb). 

3. Results and Pi scussion 


3.1 Summary of results of [1] : Key results from [1], for the "LRN" airfoil, 

are shown in Figures 2 to 5. The wake velocity spectra for different a are shown 
in Figure 2. At large a spectral peaks are observed at relatively higher fre- 
quencies which correspond to St s » 0.2. This is due to the familiar bluff-body 
shedding involving the asymmetric Karman vortex street. Around the onset of static 
stall, however, at about a = 15° in this case, the low frequency oscillation is 
observed, the corresponding St s being an order of magnitude lower. It is this 
low St s oscillation — its origin, significance as well as its contrast to the 
bluff-body shedding — that has been addressed in [1], 

The occurrence of the low St s oscillation was found to be rather illusive 
and sensitive to ambient conditions. This is delineated by the data in Figures 3 
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and 4. The lower pair of spectra in Figure 3 shows that with the normal tunnel 
operation yielding a turbulence intensity of about 0.1 percent, the periodic oscil- 
lation does not take place; however, it does when the free-stream turbulence is 
raised with the help of the screen inserted upstream of the test section. Note 
that the oscillation was observed previously in a Langley tunnel in which the free- 
stream turbulence was also higher — about 0.25 percent [2]. The upper pair of 
spectra in Figure 3 demonstrate that the low frequency oscillation can also be pre- 
cipitated if acoustic excitation in a certain high frequency range is applied. 

On the other hand, when the oscillation takes place naturally, e.g., in the 
flow with the higher free-stream turbulence, acoustic excitation in a relatively 
higher frequency range is found to eliminate it. An example is shown by the pair 
of spectra in Figure 4. The ranges of excitation frequency augmenting or suppress- 
ing the St s » 0.02 oscillation have been documented for various conditions in 
[1], The higher frequency excitation producing the suppression effect shown in 
Figure 4 appears to involve "acoustic tripping" of the separating boundary layer. 

It is conjectured that the boundary layer around the region of separation is 
turned fully turbulent by the acoustic excitation — in which case the low fre- 
quency oscillation does not take place. Thus, on the one hand the phenomenon does 
not seem to take place in clean flows yielding laminar boundary layers, while on 
the other, it also disappears when the separating boundary layer is fully turbu- 
lent. Therefore, a certain transitional state near the separation region seems to 
be a requisite for the low frequency oscillation to take place [1]. 

Also included in [1] are the computational results for the "LRN" airfoil. 

C.L. Rumsey had previously observed a similar low frequency oscillation computa- 
tionally for a NACA001 2 airfoil at a = 18° [5]. The computations were then car- 
ried out for the Reynolds number and airfoil shape as in the experiment. The 
results were found to capture the essential features of the experimental observa- 
tion. The computed flow field agreed in overall features with corresponding 
measurements. The computed force fluctuations also agreed generally with the 
experimental results, as shown in Figure 5. Despite the difference in the 
Strouhal number, the amplitudes and the details of the fluctuations are found to 
have close similarities. Aspects of Rumsey's computation will be discussed fur- 
ther i n the fol lowi ng . 

3.2 Present experimental results : For comparison with the "iced" airfoil, the 

lift, drag and moment coefficient variation with a for the "LRN" airfoil are doc- 
umented in Figure 6. Note that the low frequency oscillation takes place around 
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15° (within a narrow range of a) where the forces and the moment go through a 
large change due to stall. 

Corresponding aerodynamic data for the "iced" airfoil are shown in Figure 7. 

In comparison to the "LRN" airfoil, the stall occurred at a much lower angle of 
attack of about 7°. No hysteresis was observed around the stall angle. It is 
slightly above this a where this airfoil exhibited a similar low frequency 
osci 1 lation. 

Wake velocity spectra similar to those in Figure 2 are presented in figure 8 
for the "iced" airfoil. The oscillation is observed to occur around the a-range 
of 8° to 9°. Note that in this case the spectral peak is not as sharp as in the 

"LRN" airfoil case, but the tendency towards an oscillation around St s = 0.02 is 

unambiguous. The oscillation at three different Reynolds numbers at a = 8.5° is 
documented in Figure 9. R c = 1 .25x1 0^ corresponds to the highest speed obtainable 
i n the tunnel . 

The effect of acoustic excitation on the low frequency oscillation over the 
"iced" airfoil was briefly studied in the present experiment. An attempt was made 
to augment the oscillation as in the "LRN" airfoil case (Figure 3). However, the 

spectral peak could not be made much sharper under the acoustic excitation. 

Increasing the free-stream turbulence with the screen also did not make a signifi- 
cant difference. However, the suppression of the low frequency oscillation, as 
shown in Figure 4 for the "LRN" airfoil, could be achieved. An example is shown 
in Figure 10. 

The data in Figures 7 to 10 lead us to believe that the oscillation taking 
place around a « 8.5° with the "iced" airfoil is morphologically the same as that 
studied in [1]. The differences observed with respect to the effect of acoustic 
excitation are not surprising and reinforces the notion that the phenomenon is 
very sensitive to ambient perturbations affecting the state of the separating 
boundary layer. 

It was observed in [1] that the phenomenon seems to occur with airfoils exhib- 
iting "trailing edge stall" or "leading edge stall accompanied by a separation 
bubble", but not with ones exhibiting the "abrupt leading edge stall". The compu- 
tations and experiments suggest that the "iced" airfoil belong to the second of 
the aforementioned stall categories (leading edge stall accompanied by a separa- 
tion bubble), thus, corroborating the above observation. Note that the "LRN" air- 
foil is a borderline case between the first and the second categories of stall. 

An Wortman airfoil, on the other hand, is clearly characterized by the third cate- 
gory of stall, and thus, does not exhibit the low frequency oscillation [1]. 
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3.3 Similar observations by others : It appears that the low frequency oscil- 

lation phenomenon has gone practically unrecognized over the several decades of 
airfoil research. To our knowledge, there are only a few references in the litera- 
ture regarding similar observations. B.M. Jones (1934; see [1]) and Farren (1935; 
[6]) apparently encountered the phenomenon. Figure 11(a) reproduces the essential 
feature of Farren's data in this connection. Farren was actually developing a 
fast response balance and experimenting with oscillating airfoils. In his report, 
also included was the normal force variation with time for a R.A.F. 28 airfoil 
held fixed at a = 14°. This airfoil was previously observed by B.M. Jones to 
undergo violent fluctuations near stalling conditions. The average period of 
oscillation (Figure 11(a)) was observed by Farren to correspond to 13 chords of 
travel. This converts to St s = 0.019. 

At NASA Lewis there is a research program addressing the problem of airfoil 
leading edge ice accretion and its remedies. The program includes computation and 
experiments on the aerodynamic characteristics of such airfoil sections both 
in-house and under grants. Bragg and Khodadoust [7], under a grant, obtained wind 
tunnel data on the aerodynamic characteristics of a NACA0012 airfoil with simulated 
ice shapes. Recently, for a "glaze ice" shape, (same as in the present study), 
they also observed a low frequency flow oscillation. Figure 22 of their paper is 
reproduced as Figure 11(b). The data represent u'-spectrum measured near the 
upper edge of the shear layer, shortly downstream of the "ice horn", for a » 6°. 
(Note that "a = 4°" quoted in their paper is measured with respect to the chord of 
the basic NACA0012 airfoil, a, measured with respect to the chord as shown i.i 
Figure Kb), is about 2° larger). There is a clear peak in the spectrum, and this 
corresponds to St s » 0.016. The authors did not find such spectral peak at other 
locations over the airfoil and, therefore, commented that it should represent a 
local shedding from the ice horn and not a global bubble oscillation. However, 
their observations were preliminary, and it would appear that the same phenomenon, 
as addressed here, has been encountered by them. Their data, however, represent a 
significantly higher R c of 1.5x1 0^ , indicating that the phenomenon may not be 
characteristic of low Reynolds number airfoils only. Other relevant observations 
in the literature, including the characteristic frequencies in dynamic stall and 
stall flutter, have been discussed in [11. 

3.4 Computational result : A two-dimensional Navi er-Stokes code coupled with 

a C-mesh grid generation code was used for the computations, the details of which 
can be found in [3]. A "modified mixing length" (MML) turbulence model was used 
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in the computation. An unsteady flow with "periodic vortex shedding" was noted 
around a = 8°. The lift coefficient variation with time, computed for the same 
flow conditions of C 3 3 , is shown in Figure 12(a). Similarities can be observed 
with the Ci variation for the "LRN" airfoil (Figure 5). However, when the Bald- 
win-Lomax turbulence model was used in the present computation the oscillations 
damped out in the resultant flow field. This is shown in Figure 12(b). 

Here, let us note that in the computations by Rumsey that resulted in low fre- 
quency oscillation for the "LRN" airfoil (Figure 5), the Bal dwi n-Lomax turbulence 
model was used. Rumsey [5] also observed similar low frequency oscillations for a 
NACAOO 1 2 airfoil at a = 18° and R c = 10 6 . 

The algorithms used in [33 and [53 are basically similar. The "MML" is an 
algebraic turbulence model which attempts to avoid some of the difficulties encoun- 
tered with Baldwin-Lomax turbulence model when applied to separated flows. 

Details of the model will be discussed in the second author's dissertation, which 
is forthcoming. The "turbulent viscosity" (pt> obtained from the two turbulence 
models are quite different. Typically, the "MML" model obtains pt^l values on 
the order of 10^ while in a comparable flow the Baldwin-Lomax model obtains values 
on the order of 10^, where y] is the laminar viscosity. It is probable that the 
"transitional state" of the separating boundary layer, conjectured in [13 to be 
necessary to produce the low frequency oscillation, may be characterized by certain 
values of the ratio y^Hl • However, a systematic trend is far from clear at this 
point and will require further investigation. 

The computations were also carried out for the "iced" airfoil, at a = 9°, 

R c = 1C>5 and M = 0.10, to approximately match the present experimental condi- 
tions. The results have added further to the confusion. The "MML" turbulence 
model, which yielded the oscillation at the higher R c , resulted in damped oscilla- 
tions. The corresponding C] variation is shown in Figure 13(b). On the other 
hand, when the flow was assumed laminar for this case, the computation yielded a 
low frequency oscillation! This is shown in Figure 13(a). The St s values corre- 
sponding to the oscillations in Figures 12(a) and 13(a) turn out to be about 0.011 
and 0.008, respectively. Let us note, for comparison, that Rumsey' s computation 
for both airfoils ("LRN" and NACA0012) yielded an St s ~ 0.03. The experimental 
and the computational results are summarized in Tables 1 and 2, respectively. 
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4. Concluding Remarks 


The experimental evidence gathered in this paper suggest that the low fre- 
quency oscillation may indeed be a frequent occurrence with various airfoils. Its 
occurrence with "glaze ice accretion" is certainly a possibility. Thus, the lead- 
ing edge ice accretion may not only be detrimental for aerodynamic performance but 
should also be of concern in terms of large unsteady loads associated with the phe- 
nomenon. The phenomenon, however, occurs in a narrow range of angle of attack 
around the onset of static stall. The acoustic excitation effect results also sug- 
gest a possible remedy, which needs to be explored further. The accumulated exper- 
imental evidence points towards a Strouhal number of about 0.02 associated with 
the phenomenon, although there are scatter in the data. With respect to the compu- 
tations, questions remain in the application of turbulence models to separated 
flows. Nevertheless, it is felt that the essence of the phenomenon can be captured 
computationally with certain combinations of turbulence model, Reynolds number and 
airfoil shape. Obviously, a detailed computational study is in order not only to 
explore the validity of the algorithms but also to shed light into the still unre- 
solved mechanism of this curious phenomenon. 
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TABLE 1 . - EXPERIMENTAL RESULTS ON THE LOW FREQUENCY 

OSCILLATION 


Ai rfoi 1 

a 

R c xl0 -3 

c, 

cm 

St s 

Refer- 

ence 

R.A.F.28 

14 

100 

15.2 

0.019 

[6] 

"LRN" 

15 

40 to 140 

10.2 

0.02 

C 2 3 

"LRN" 

17 

15 to 25 

7.3 

0.025 

[1] 

"LRN" 

15 

50 to 150 

12.7 

0.02 

[1] 

"LRN" 

16 

100 to 300 

25.4 

0.033' 

Cl] 

NACA001 2 

14 

70 

10.2 

0.022 

[1] 

(w/trip) 

Wortman 

15 

50 to 150 

12.7 

Damped 

[1] 

"Iced" 

6 

1500 

55 

0.016 

[7] 

"Iced" 

8.5 

75 to 125 

12.7 

0.02 

Present 


fHigher mean velocity due to large blockage may have 
resulted in higher St s . 


TABLE 2. - COMPUTATIONAL RESULTS 


Ai rfoi 1 

a 

Turbulence 

model 

R c x 1 0~ 3 

M 

St s 

Reference 

NACA001 2 

18 

Baldwin-Lomax 

1000 

0.3 

0.03 

[5] 

NACA001 2 

18 

Lami nar 

1000 

0.3 

0.16 

C 5 3 

"LRN" 

15 

Baldwin-Lomax 

75 

0.3 


[1] 

Wortmann 

mm 

Baldwin-Lomax 

75 

0.3 

Damped 

Cl] 

"Iced" 

ipi 

Baldwin-Lomax 

1400 

0.12 

Damped 

Present 

"Iced" 

H 

MML 

1400 

| 

0.12 

0.011 

Present 
and [3] 
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MML 

100 

0.10 

Damped 

Present 


9 

Lami nar 

100 

0.10 

0.008 

Present 
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FIGURE la. - WIND TUNNEL TEST SECTION. 


FIGURE 1b. - AIRFOIL CROSS-SECTIONS. UPPER: 
LRN(1)-1007, LOWER: NACA0012 WITH 'GLAZE 

ICE ACCRETION'. 



FIGURE 2. - WAKE VELOCITY SPECTRA FOR 
'LRN' AIRFOIL. R c = 10 5 , u'^/U^ 
O.A PERCENT; MEASUREMENT AT x/C = 1.5, 
y/C = 0.15, Z = 0. SPECTRA TRACES 
ARE STAGGERED SUCCESSIVELY BY ONE OR- 
DINATE DIVISION. 
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FIGURE 3. - WAKE VELOCITY SPECTRA FOR 
'LRN' AIRFOIL AS IN FIGURE 2, SHOW- 
ING ENHANCEMENT OF St s * 0.02 OSCIL- 
LATION. R c = 10 5 , a = 15°. UPPER 
PAIR OF TRACES SHIFTED BY TWO ORDI- 
NATE DIVISIONS. 



FIGURE 4. - WAKE VELOCITY SPECTRA FOR 
'LRN' AIRFOIL AS IN FIGURE 2 . SHOW- 
ING SUPPRESSION OF St $ « 0.02 OSCIL- 
LATION. R c = 0.75 x 10 5 , a = 15°. 
u'oo/Uoo= 0.4 PERCENT. 




FIGURE 5. - LIFT-COEFFICIENT VARIATION 
WITH TIME FOR 'LRN' AIRFOIL AT a = 15° 
AND R c = 0.75 x 10 5 . BALDWIN-LOMAX 
TURBULENCE MODEL USED IN COMPUTATION. 
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FIGURE 8.- WAKE VELOCITY SPECTRA AS IN 
FIGURE 2 FOR THE 'ICED' AIRFOIL; R c 
= 10 5 . TRACES ARE STAGGERED BY ONE 
DIVISION. 



FIGURE 9. - WAKE VELOCITY SPECTRA AS 
IN FIGURE 2 FOR THE 'ICED' AIRFOIL 
AT a= 8.5° FOR INDICATED R c . TRA- 
CES ARE STAGGERED BY ONE MAJOR DIV- 
ISION. 
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FIGURE 10.- WAKE VELOCITY SPECTRA AS IN 
FIGURE 4 FOR THE 'ICED' AIRFOIL. SHOW- 
ING SUPPRESSION OF St s * 0.02 OSCIL- 
LATION. a = 8.5°. R c = 10 5 . 


FIGURE 11. - (A) NORMAL FORCE VAR- 
IATION FOR R.A.F. 28 AIRFOIL AT 
a = 14° AND R c = 10 5 FROM REF. 161. 
(B) VELOCITY SPECTRA MEASURED SHORT- 
LY DOWNSTREAM OF THE 'HORN' OF AN 
' ICED' AIRFOIL AT a ® 6° AND R c = 
1.5 X 10 6 FROM REF. 171. 
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(b) BALDW IN-LOMAX TURBULENCE MODEL. 



(a) "MML" TURBULENCE MODEL. 

FIGURE 12. - COMPUTED LIFT-COEFFICIENT 
VARIATIONS WITH TIME FOR THE 'ICED' 
AIRFOIL AT A HIGHER R c . R c = 1.M X 
10 6 , M = 0. 12/ a = 10°. 




TUoo/c 

(a) LAMINAR FLOW. 


FIGURE 13. - COMPUTED LIFT-COEFFICIENT 
VARIATIONS WITH TIME FOR THE 'ICED' 
AIRFOIL AT THE R c COVERED IN THE EX- 
PERIMENT. R c = 10 5 , M = 0.1. a = 9°. 
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